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Motivated by experiments in which a polynucleotide is 
driven through a proteinaceous pore by an electric field, we 
study the diffusive motion of a polymer threaded through 
a narrow channel with which it may have strong interac- 
tions. We show that there is a range of polymer lengths in 
which the system is approximately translationally invariant, 
and we develop a coarse-grained description of this regime. 
From this description, general features of the distribution of 
times for the polymer to pass through the pore may be de- 
duced. We also introduce a more microscopic model. This 
model provides a physically reasonable scenario in which, as 
in experiments, the polymer's speed depends sensitively on 
its chemical composition, and even on its orientation in the 
channel. Finally, we point out that the experimental distribu- 
tion of times for the polymer to pass through the pore is much 
broader than expected from simple estimates, and speculate 
on why this might be. 



INTRODUCTION 

Modern polymer physics has achieved great success 
with models in which the polymer is regarded as a flexi- 
ble, uniform "string" whose conformational entropy dom- 
inates the system's behavior (de Gennes, 1979; Doi and 
Edwards, 1986). Although this is usually an excellent 
description, in some situations other interactions can be- 
come important. One example is the insertion of a poly- 
mer into a pore of diameter comparable to the size of 
the chemical repeat units that make up the polymer. Al- 
though perhaps unusual with synthetic polymers, such 
a situation can easily occur in biological systems. For 
example, Kasianowicz, Brandin, Branton, and Deamer 
(hereafter KBBD) have recently detected single strands 
of RNA (polyuridylic acid) passing through a 1.5 nni 
pore formed by a membrane-bound protein (Kasianow- 
icz et al., 1996). Szabo and coworkers (Szabo et al., 
1997; Szabo et al., 1998) and Hanss and coworkers (Hanss 
et al., 1998) have studied similar systems. In addition to 
their intrinsic interest, these experiments may eventu- 
ally lead to a single-molecule RNA and DNA sequenc- 
ing technique. More generally, most cells must transport 
macromolecules across membranes in order to function; 
in several cases, relatively "thick" molecules are believed 
to pass through nanometer-scale channels. The translo- 
cation of polynucleotides through proteic pores has been 
implicated in a variety of processes, including phage in- 
fection and bacterial conjugation (Dreiseikelmann, 1994), 
the uptake of oligonucleotides by certain organs (Hanss 
et al., 1998), and transport across the nuclear envelope 



in plants (Citovsky and Zambryski, 1993). It has been 
speculated that some of these transport pathways could 
eventually prove important in gene therapy (Szabo et 
al., 1998; Hanss et al., 1998). Similarly, polypeptide- 
conducting channels play an important role in protein 
kinesis (Schatz and Dobberstein, 1996; Simon and Blo- 
bel, 1991); in a few instances, the translocation may even 
be driven by electrophoretic effects (Attardi and Schatz, 
1988). 

There exists a considerable literature on the confine- 
ment of polymers in channels of diameter significantly 
larger than the polymers' persistence length (de Gennes, 
1979); well-developed scaling techniques can be used in 
the theoretical treatment of this regime. Recently, theo- 
rists have also shown an interest in the opposite limit of a 
very narrow, almost point-like hole. For example, several 
groups have studied the diffusion of polymers across ide- 
alized, infinitely thin membranes (Carl, 1998; Di Marzio 
and Mandell, 1997; Yoon and Deutsch, 1995; Lee and 
Obukhov, 1996; Park and Sung, 1998a; Park and Sung, 
1998b; Sung and Park, 1996). The pore and the mem- 
brane are viewed as hard walls whose only interaction 
with the polymer is steric, and the emphasis is on how the 
walls' presence decreases the entropy and slows the dy- 
namics of those parts of the polymer outside of the hole. 
Possible mechanisms for the active transport of polymers 
through pores in biological systems have also been stud- 
ied (Peskin et al., 1993; Simon et al., 1992; Sung and 
Park, 1996). 

Inspired largely by the experiments of KBBD, in this 
paper we consider an different scenario: We study the 
motion of a homopolymer threaded through a narrow 
pore with which it has strong interactions. The pore 
is taken to be sufficiently small that no more than one 
polymer diameter can fit in it at a given time; in par- 
ticular, "hairpin" bends are not allowed to pass through 
the channel. We also put aside the question of how the 
polymer first enters the hole, focusing instead on the dy- 
namics once one end has been inserted. We then ar- 
gue that, in the presence of a force driving the polymer 
through the pore, there should be a regime in which the 
polymeric degrees of freedom outside of the pore can be 
neglected, and the system is effectively one-dimensional. 
In this limiting case we propose a two-tiered picture: a 
coarse-grained "macroscopic" description of wide valid- 
ity and a simple "microscopic" model from which the 
"macroscopic" parameters may be calculated. Our ap- 
proach follows several authors (Peskin et al., 1993; Si- 
mon et al., 1992; Park and Sung, 1998a; Park and Sung, 
1998b; Sung and Park, 1996) in viewing the transloca- 



tion process as essentially diffusion in one dimension; we 
differ, however, in emphasizing the role that interactions 
with the pore itself play in this diffusion process. On the 
more microscopic level, we include the effects of these 
interactions through a tilted washboard potential, simi- 
lar to models of laser mode locking (Haken et al., 1967) 
or phase dynamics in Josephson junctions (Ambegokar 
and Halperin, 1969) (see figures g and ^). The periodic 
modulation of the potential refiects the periodicity of the 
polynucleotide's sugar-phosphate backbone. The impor- 
tance of polymer-pore interaction has previously be em- 
phasized by Bezrukov and coworkers (Bezrukov et al., 
1996); our model also bears some similarity to work on 
gel electrophoresis that examines the importance of local 
"solid friction" forces between the polyelectrolyte and the 
gel (Deutsch, 1987; Burlatsky and Deutch, 1993; Viovy 
and Duke, 1994; Burlatsky and Deutch, 1995; Deutsch 
and Yoon, 1997). Although the macroscopic parame- 
ter values for KBBD's system differ in some respects 
from those predicted by our microscopic model, we are 
nonetheless able to make several fairly robust predictions. 
More importantly, we show how a simple physical mecha- 
nism can account for several striking features of the data 
of KBBD. We thus hope that our work will provide a 
useful contribution to our understanding of the translo- 
cation of polyelectrolytes. 

Since our analysis relies heavily on KBBD's results, 
the next section sketches some salient features of their 
data. We then introduce a long length scale "hydrody- 
namic" description of one-dimensional diffusion and use 
it to calculate the distribution of passage times for a poly- 
mer being driven through a pore. The arguments used 
to arrive at these results are quite general; in particular, 
they require few assumptions about the details of the mi- 
croscopic dynamics of the system. There are, however, 
circumstances when our approximations break down, and 
we consider these next. This section also serves to em- 
phasize several aspects of the experiments that will guide 
our choice of the microscopic model in the succeeding sec- 
tion. After introducing this microscopic model, we use 
it to calculate a mean drift velocity and an effective dif- 
fusion coefficient and compare them to values estimated 
from KBBD's data. These comparisons will reveal cer- 
tain features that cannot be accounted for by our model 
in its simplest form, so we then discuss possible reasons 
for the discrepancy, as well as touching brieffy on several 
applications of our calculations. We conclude by sum- 
marizing our results and highlighting some issues that 
remain open. 



EXPERIMENTAL BACKGROUND 

In the experiments of interest to us, KBBD worked 
with a Staphylococcus aureus a-hemolysin ion channel in 
an artificial lipid bilayer membrane (diphytanoyl-PC). 
This channel has the advantage that for many pur- 



poses it may be considered always to be open. The a- 
hemolysin protein has recently been crystallized and an 
x-ray structure obtained (Song et al., 1996). This reveals 
a mushroom-shaped complex with a roughly 10 nm long 
solvent-filled channel. The channel is 1.5 nm in diameter 
at its narrowest constriction, barely larger than the diam- 
eter of a single polynucleotide strand. After inserting a 
single pore into a bilayer membrane and applying a trans- 
membrane potential of between 110 and 140 mV, KBBD 
added homopolymeric single- stranded DNA or RNA to 
one side of the membrane, designated cis. The samples 
of polynucleotides had mean lengths on the order of a few 
hundred nucleotidea] and were assumed to be close to 
monodisperse. After adding the polynucleotides, KBBD 
monitored the transmembrane ionic current as a function 
of time. The time series shows a baseline current, modu- 
lated by periods on the order of hundreds of microseconds 
in which the current decreases almost to zero (figure n, 
inset). A variety of observations support the interpreta- 
tion that these blockades were caused by the passage of a 
polymer through the a-hemolysin channel. The data of 
KBBD can thus be interpreted as giving measurements 
of the times required for individual polynucleotides to 
traverse the membrane under the infiuence of an electric 
field. 

When these data are displayed as a histogram, with 
the number of observed events plotted against the length 
of the blockade (figure ||), one sees that the blockade 
times fall into three distinct peaks. Of these, the first 
("peak 1") is caused by polymers that enter and retract 
and thus do not completely cross the membrane, while 
the other two ("peak 2" and "peak 3") are both the re- 
sult of a polymer's actually passing through the channel. 
The polymers in peak 3 evidently cross the membrane 
roughly three times faster than those in peak 2. KBBD 
made the intriguing suggestion that there are two charac- 
teristic times associated with translocation because the 
polynucleotide can enter the pore in two distinct direc- 
tions: One peak corresponds to polymers that enter the 
channel with their 3' end first, the other to polymers that 
enter with their 5' end first. We will show in subsequent 
sections how such behavior can arise from a simple mi- 
croscopic model. 

A quantity of considerable interest in what follows will 
be the mean force F driving the polymer through the 
porcDJ Clearly F is primarily the result of the electric 



'Various groups have measured the persistence length of 
single-stranded DNA in high salt concentrations to be be- 
tween 0.75 nm and 1.5 nm (Achter and Felsenfeld, 1971; Smith 
et al., 1996; Tinland et al., 1997), or roughly 1 to 2 nu- 
cleotides, meaning that the polymers used were of order 100 
persistence lengths long. 

^ One can define F more precisely as the mean force required 
to immobilize a given monomer in the pore, where the av- 
erage is taken over time and over all of the monomers in 
a given polymer. Thus F does not include hydrodynamic 



field acting on the polymer. Since a long, narrow channel 
has a much larger electrical resistance than the macro- 
scopic volumes of solution on either side of the mem- 
brane, any voltage V applied to the system should fall al- 
most entirely across the a-hemolysin pore. The charge on 
each nucleotide is just the electron charge e, so the elec- 
trostatic energy gained by moving one nucleotide through 
the pore is eV. This suggests that F is roughly 



a 



.knT 



(1) 



where a « 6 A is the length of a nucleotide, and the second 
equality holds for V ~ 125mV. For most of the rest of 
the paper, we will assume that F takes this value, and 
many of our arguments will be based on the fact that F 
is thus quite large when expressed in appropriate units. 
Some effects that could modify F are considered in the 
discussion section and in appendix C. 

Before presenting our model, we would finally like to 
review the experimental evidence that the interactions 
between the polymer and the a-hemolysin pore do in- 
deed play the dominant role in KBBD's experiments. 
We have already mentioned the existence of two dis- 
tinct characteristic times for the polymer to cross the 
membrane. Such a result is easiest to interpret if one 
believes that the polymer's speed is determined by in- 
teractions between the polymer and the narrow channel 
constriction, where molecular scale asymmetries could 
be important. Similarly, recent data show that homo- 
polynucleotides of different bases can move at strikingly 
different speeds (Dan Branton, Harvard University, per- 
sonal communication): poly[U] is of order 20 times faster 
than poly[dA]. Although chemical differences certainly 
can lead to variations in polymer properties such as the 
persistence length, we believe that such strong depen- 
dence on molecular details can more easily be explained 
if we focus on the pore region. Finally, even the fastest 
polynucleotides pass through the pore far more slowly 
than simple estimates of hydrodynamic drag would sug- 
gest: Model the pore as a cylindrical hole of radius R 
and the part of the polymer in the pore as a cylinder 
of radius r. Then, when the polymer moves with speed 
V, the drag force per length on the part in the pore is 
roughly 2Tri]rv/{R — r). Electrophoretic effects change 
this result very little (see appendix C). For a polynu- 
cleotide in an a-hemolysin channel, r/(R — r) is some- 
what larger than unity, and the total length of the "cylin- 
der" is roughly 50A. According to scahng arguments of 
Lee and Obukhov (Lee and Obukhov, 1996), the contri- 



drag forces, nor forces that vanish when averaged over all the 
monomers. Equivalently, F can be defined by requiring that 
exp{Fa/kBT) be the ratio of the probabilities that the poly- 
mer will move forward one base to the probability that it will 
move backwards one base, again appropriately averaged over 
all monomers. 



bution to the drag force from the ends of the polymer 
outside the channel is only 2 x 6TTr]bv, where rj is the 
solvent viscosity, and the Kuhn length b is between 15 
and 30A. Even if hydrodynamic interactions are entirely 
screened by the motion of counterions (as they are for the 
electrophoresis of an isolated polymer in solution, with 
screening length of order the monomer size) , the drag on 
those parts of the polymer in solution cannot be larger 
than roughly Air-qLv. If one substitutes typical parame- 
ter values for KBBD's experiments and balances the sum 
of these drag forces with the naive driving force of ^k-gT 
per nucleotide, one finds that the polymer would be ex- 
pected to move through the pore at a rate of roughly 
lO^nucleotides/second, 100 times faster than observed. 
The three observations of this paragraph, taken together, 
certainly suggest that we focus on the degrees of freedom 
in the pore when trying to understand the experiments 
of KBBD. 



COARSE-GRAINED DESCRIPTION 

Motivation and Governing Equation 

This section, and most of the rest of the paper, is con- 
cerned with predicting distributions of blockage times of 
the sort shown in figure H. It is now well-established in 
condensed matter physics that the form of the slow, long 
length-scale dynamics of a system is often determined by 
the system's symmetries and conservation laws. All mi- 
croscopic details are subsumed in phenomenological cou- 
pling constants and transport coefficients. In this spirit, 
we would like to obtain a coarse-grained equation for the 
probability P{x,t) that a contour length x of the poly- 
mer's backbone has passed through the pore at time t. 
(The variable x is defined so that if the polymer backbone 
has length L, x = when the polymer has just started in 
the pore and x = L when it has reached the other side). 
For such a "hydrodynamic" description to make sense, 
several conditions must be met. One is that the polymer 
length L be much larger than the distance a between suc- 
cessive nucleotides. We also demand that the dissolved 
counterions (as well as the solvent and any other solutes 
) relax quickly compared to the translocating polymer, 
in order that we may ignore their dynamics. Since the 
ions are much smaller than a polynucleotide, and conse- 
quently diffuse much faster, this condition should not be 
difficult to satisfy. Finally, our task will be considerably 
simplified if the microscopic system is (approximately) 
invariant under translations by an integer multiple of a in 
either direction. Then, after averaging over variations on 
the scale of a single nucleotide, we must obtain a transla- 
tionally invariant equation. We will give this assumption 
a firmer basis in the next section. Roughly, however, 
there should be translational invariance when we can ne- 
glect the parts of the polymer outside of the channel, 
and this in turn should be possible when the interactions 



between the polymer and the pore are strong enough. 

Under the conditions just outhned, the (probabihty) 
density of the polymer is the only conserved variable, and 
it is relatively straight-forward to write down the coarse- 
grained hydrodynamic equation for P. Because there is 
only a single polymer (or, equivalently, a "gas" of non- 
interacting polymers going through the same hole), the 
probability current j, defined by dP/dt+dj /dx = must 
be linear in P. The lowest order allowed terms are then 
proportional to P and to dP/dx: 



j{x,t)=vP{x,t)-D 



dP{x,t) 
dx 



(2) 



The first term is permitted because there is an electric 
field driving the system. P then satisfies the familiar 
equation for diffusion with drift, 



dP _ d^P 
dt dx'^ 



dP 



(3) 



Here v and D are respectively an average drift velocity 
and an effective diffusion coefficient. Their values are 
determined by more microscopic physics; in particular, 
they may depend nonlinear ly on the applied electric field. 
Eq. H may alternatively be derived from a microscopic 
master equation that is invariant under translations by 
a. The coefficients v and D are then related to the lowest- 
lying eigenvalues of the master equation. This connection 
will be illustrated in a subsequent section. 

On the macroscopic level of Eq. H, all information on 
the competition between driving and diffusive spreading 
is encoded in a parameter that we call the diffusive length 
1,1 = D/v. Roughly speaking, on length scales less than 
Id^ the polymer's motion is little affected by the pres- 
ence of the bias from the electric field, while on scales 
larger than Id, the driving dominates. Indeed, if Eq. 
described a rigid particle diffusing in one dimension un- 
der the infiuence of a uniform force /, an Einstein re- 
lation would hold, and we would have v = Df /(UbT), 
and Id — k^T / f . Thus, in this case. Id is precisely the 
length over which the driving force does a quantity k^T 
of work. In the remainder of this section, we will often 
assume that the length L of the polymer is larger than 
Id, a condition satisfied by KBBD's data. 



Distribution of Passage Times 

We now propose to calculate a distribution of passage 
times of the sort measured by KBBD. This section will 
show that, for given v and D, the probability ■ip{t) that 
the polynucleotide takes a time t to pass through the 
channel has only one peak. Thus, the presence of two 
peaks in KBBD's data must be explained by the assump- 
tion that different physical situations give rise to different 
values of v and D. Subsequent sections will argue that a 
polynucleotide passing through the pore with its 3' end 



first can indeed have an average velocity that is signifi- 
cantly different from one passing through with its 5' end 
first. This section, however, is confined to the calcula- 
tion of the passage times for fixed parameter values. The 
distribution ^{t) we obtain should thus be compared to 
a single peak in the data of KBBD. 

One can easily estimate the first few cumulants of this 
distribution. If a polymer of length L moves with aver- 
age velocity u, one expects that the mean time to pass 
through the channel should be (t) ~ L/v. Likewise, the 
variance in the distance traveled in a time {t) is (Ax)^ = 
2D{t). It would then seem reasonable that the variance 
in arrival times should be At^ = ((i- (i))^) « (Ax)^/w^, 
or At^ « 2DL/v^. These conclusions are in fact roughly 
correct for a sufficiently long polymer. One might expect 
corrections, however, because some fraction of the poly- 
mers that enter the pore will leave again from the same 
side instead of passing all the way through. On average, 
these will be the "slower" molecules: Those that spend a 
significant time with only the tip of the polymer inserted 
in the channel are far more likely to fall back out than 
are those that are quickly driven through the hole. Thus, 
only "faster" chains tend to enter into the calculation of 
the mean transit time, decreasing (i). This effect is most 
pronounced for small L/ld, because only molecules within 
Id of the cis side have an appreciable chance of "backing 
out" instead of exiting on the trans side. Indeed, when 
L <ti Id, the driving should be negligible, and we expect 
(t) to approach its u = value L^/QD. To determine the 
precise form of this crossover, we must turn to a more 
detailed calculation. 

This calculation can be formulated as one of a well- 
studied class of problems known as first-passage prob- 
lems (Risken, 1984; van Kampen, 1992). Essentially, all 
that is required is to solve Eq. 13 on the interval [0, L] 
with absorbing boundary conditions P(0) = P{L) = 0. 
Then, the current density j{L) at L gives the probability 
per time that the polymer will leave the pore from the 
far (trans) side, while — j(0) is the probability per time 
that it will exit from the cis side from which it entered. 
One must also specify the starting point xo G [0, L] of 
the polymer; in what follows, we always take the limit 
xq — *■ 0, in keeping with the fact that the polymer starts 
entirely on the cis side of the membrane. The algebraic 
details of the solution are summarized in appendix A; 
here we include only a discussion of the main results. 

Although exact expressions for (t) and At may be ob- 
tained, it turns out to be more instructive to consider 
the distribution ip{t) itself. For arbitrary L/ld, this can 
only be expressed as an infinite series, but if terms that 
become exponentially small as L"^ /(vtld) — > cxd are ne- 
glected, a comparatively simple analytic expression is ob- 
tained: 




idivtfi^ [vtyi^ 



^-(vt-Lf/{ivtU) 



(4) 



Note that this expression is not vahd for sufficiently large 
t, and in particular not for t so large that it predicts that 
^p{t) becomes negative. Nonetheless, for values of t near 
the maximum in V'(^)i i-S- those such that vt/L ~ 0(1), 
it is accurate to within a percent for L/ld as small as 
4, and correctly reflects the qualitative features of ip{t) 
for significantly smaller L/ld- Figure g plots ip(t) for 
L/ld = 5; a Gaussian with the same mean and vari- 
ance is included for comparison. Evidently, ipit) is quite 
skewed, and its mean and maximum are correspondingly 
well-separated. Thus, (t) and At are not the best pa- 
rameters for describing experimental data. Indeed, both 
cumulants are sensitive to how 'ip{t) decays for large i, 
making them very hard to extract accurately from real- 
istic data sets. A more useful choice of parameters to 
characterize ip{t) are the position imax of its maximum 
(which satisfies dip/dt\f = 0) and the width St of the 
peak. The latter is defined as St = (tn — iL)/2, where tu 
and II satisfy ^'(i-Ri^i) = e~^'^'0(imax); we have chosen 
a factor of e~^/^ instead of the more conventional 1/2 to 
facilitate comparison with fits of data to a Gaussian. One 
expects that as L/ld — * oo, tmax and St should approach 
(t) and At, respectively. For example, for large L/ld we 
have, 
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(5) 



The rapidly growing coefficients indicate that although 
^max approaches L/v as L approaches infinity, it falls 
away from its asymptotic form quite rapidly for finite 
L. 

More generally, one can easily find tmax and St by nu- 
merically solving the equations that define them. Fig- 
ure plots (5t/tinax versus the polymer length L. This 
ratio is especially interesting because it depends only on 
L/ld^ and not on v and D separately; one can thus use 
it quickly to estimate L/ld- In KBBD's data, St/tj^^x is 
usually of order 0.5 for a ~ 200 nucleotide chain, suggest- 
ing that L/ld ~ 5, or that Id is of order 40 nucleotides. As 
figure indicates, in this range tmax already deviates sig- 
nificantly from the naive guess tmax ~ L/v. In particular, 
tmax/L varies by a factor of 2 as L/ld increases from 5 
to 25. With sufficiently good data, this deviation from a 
strict proportionality to L might well be observable, pro- 
viding strong confirmation of our quasi one-dimensional 
picture. 



REGIME OF VALIDITY 



piece that is actually inside the channel, and the two 
"ends" , comprising the majority of the nucleotides, out- 
side the channel. The pore always contains the same 
number of bases, so, for the homopolymers, this part of 
the polymer always satisfies the requirement of trans- 
lational symmetry. The length of each end "dangling" 
outside the pore, in contrast, changes with the transloca- 
tion parameter x, destroying translational invariance. In 
what follows, we shall argue that under certain conditions 
this variation may be neglected. Our arguments assume 
that the parts of the polynucleotide outside the pore may 
be described by the theories usually applied to long, flex- 
ible polymers (de Gennes, 1979; Doi and Edwards, 1986); 
we thus ignore, for example, hydrogen-bonding and other 
specific interactions (Cantor and Schimmel, 1980). We 
also assume that the ion channel is sufficiently long and 
narrow that any voltage drop falls almost entirely across 
the channel (see appendix C). The electric field and the 
solvent flow velocity outside of the channel can then be 
ignored. 

There are two criteria for ignoring the ends of the poly- 
mer outside of the pore. First, they should have a char- 
acteristic relaxation time that is much faster than the 
characteristic time for the motion of a monomer through 
the channel. In the absence of interactions between the 
polymer and the pore, one would expect diffusion on the 
scale of a few monomers to be much faster than the re- 
laxation of a long polymer coil, and this inequality could 
never be satisfied. However, since the nucleotides in the 
pore can be expected to interact strongly with the con- 
fining protein, the requirement is not implausible. The 
longest time scale of an isolated polymer in solution is 
the Zimm time tz « OAr]R%/{kuT) « OAr^N^^b^/iltBT), 
where z/ is the Flory exponent, El b is the Kuhn segment 
length (equal to twice the persistence length), 77 is the 
solvent viscosity, and N = L/b. Substituting in numeri- 
cal values for a single-stranded polynucleotide in water, 
one finds that tz ~ A^^'^(3.2 x 10^*/isec). If we imagine 
that the polymer moves a monomer through the chan- 
nel by hopping over an energetic barrier (an idea to be 
considered in more detail when we introduce our micro- 
scopic model), then in the limit of strong driving, the 
translocation speed is simply v = a/tporc, where tporc is 
the longest relaxation time of the part of the polymer 
in the pore. Substituting numerical values for poly[U], 
we find tpore = a,/v ~ 1.5/isec. Comparing this figure to 
tz, we see that the two become of the same order when 
N is of order 150, corresponding to a length of polymer 
of roughly 300 nucleotides protruding from each side of 
the pore. Of course, for polymers that traverse the mem- 



In the previous section, we argued that a requirement 
for the validity of a one-dimensional diffusion model was 
that the system be (approximately) unchanged if the 
polymer moves an integer number of monomers forwards 
or backwards in the pore. This section discusses when 
this condition is satisfied. We begin by dividing the poly- 
mer into three parts: the roughly ten nucleotide long 



*In principle, v ~ 0.6 for a long polymer in a good sol- 
vent. However, even with the longest available chains, f is 
never observed experimentally to be larger than 0.55 (Doi 
and Edwards, 1986) , so we use this value for specific numeri- 
cal calculations. 



brane more slowly, as is the case for poly[dA], the value 
of N above tz ^ ipore can be significantly larger. 

As long as the dynamics of the polymer outside of 
the pore are fast compared to the dynamics in the pore, 
one need not treat the external degrees of freedom ex- 
plicitly. Instead, they affect the motion of the polymer 
only through a contribution J-'{x) to its free energy and 
through the increased drag they contribute.ll Lee and 
Obukhov's scaling argument (Lee and Obukhov, 1996) 
implies that their effect on the drag is independent of 
the length of polymer on a given side of the membrane. 
On the other hand, in order for us to be able to ne- 
glect JF, dT /dx must be small compared to the force F 
driving translocation. Denote the free energy of the coil 
on the cis side of the membrane by J^c{x) and that of 
the coil on the trans side by !Ft{x)] their sum is ^-"(2;). 
Sung and Park pointed out that Tc and Tt are simply 
the free energies of a polymer grafted by one end to a 
planar surface (Sung and Park, 1996). For a polymer 
of length X, this entropic free energy is known to be pro- 
portional to fcBT'ln(x/6), with a coefficient of order unity 
that depends on whether excluded volume effects are im- 
portant (Binder, 1983). Ignoring the few monomers actu- 
ally in the channel, the lengths of polymer on the cis and 
trans sides of the barrier are x and L — x, respectively, 



T{x) oc fceT 
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L — x 



(6) 



For a chain that is a fixed fraction of the way through 
the hole (i.e. for fixed x/L), dT jdx vanishes like 1/L. 
Further, it makes little sense to consider x < o, where a 
is the length of a single monomer, so we must always have 
dJ- jdx ^ k^T /a. Typical values will be much smaller 
than this bound. The driving force F « ^k^T / a thus 
greatly exceeds dT ldx\ indeed, since the polymers used 
by KBBD are several hundred nucleotides long, F is more 
than a factor of 100 larger than a typical value of dT jdx. 
In sum, we have shown that in the window of polymer 
lengths 
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the polymer is short enough to relax quickly, but long 
enough that the entropic barrier to crossing the mem- 
brane is not too steep. For lengths in this window, the 
ends of the chain hanging outside of the pore can be ne- 
glected compared to the monomers inside the pore. Since 
the system studied by KBBD falls within this window. 



we are justified in using simple one-dimensional models 
to describe it. 



MICROSCOPIC MODEL OF THE PORE 

Until now, we have avoided specifying the physics of 
the interactions within the pore. In this section, we 
present a simple phenomenological model of these in- 
teractions. Our main goal is to understand physically 
how the parameters v and D can vary sufficiently to ex- 
plain experimental facts like the difference in velocities 
between polymers moving forwards and backwards. 



Description of the Model 

We begin by focusing on the polymer backbone, whose 
coordinate x tells us what fraction of the polymer chain 
has passed through the channel. If the motion of the 
backbone is sufficiently slow compared to all the other 
degrees of freedom in the pore, then we can take x to 
be the only dynamical variable in the problem. The re- 
maining degrees of freedom are then described by a free 
energy <I>(x) that depends on the polymer translocation 
parameter x. The potential ^(x) can, for example, be 
expected to have contributions from electrostatic inter- 
actions between the polymer and the a-hemolysin hep- 
tamer. Two unit charges separated by lA in water have 
an energy of about 6kBT at room temperature; since 
both polynucleotide and protein have completely ionized 
groups in physiological pH, it is thus plausible that typi- 
cal values of $ should be at least on the order of several 
fceT. We split <f> into a mean slope F determined by 
the applied voltage drop and a part U{x) that captures 
the details of the pplvmer's interactions with the pore: 
^(x) = U{x) — Fa;.H For homopolymers (provided we 
continue to neglect the degrees of freedom outside the 
pore), U{x) is periodic, with period a = 1 nucleotide. F 
is precisely the mean force introduced in Eq. m, and is 
equal to eV/a in the simplest picture. 

Our problem is now formally no different from that of 
a point particle diffusing in a periodic potential U and 
driven by a constant force F. The probability P{x) of 
finding such a particle at a point x is governed by a 
Smoluchowski equation, 
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^Here, we assume that v is sufficiently small that the parts 
of the polymer outside the pore are essentially in equilibrium. 
On purely dimensional grounds, this must be true when tz ^ 
N^{b/v) for some non-negative exponent y, a requirement 
that is met in KBBD's experiments. 



**In principle, U could depend on the applied voltage (and 
hence on F). We ignore this effect; many of our conclusions 
will in any case turn out to be insensitive to it. 



The "bare" diffusion constant Dq is related through an 
Einstein relation to some suitable drag force on the poly- 
mer in the channel. It is not to be confused with the 
"effective" diffusion constant D that includes the effects 
of U and describes the polymer's motion on length scales 
much larger than a. 

It is helpful both for numerical work and for intuition 
building to have a concrete idea of the simplest form U{x) 
could take. In particular, such a simplified "cartoon" will 
give us an idea of the minimum number of parameters 
needed to describe the gross features of the potential. A 
natural choice for such a U{x) is a sawtooth potential of 
the sort sketched in figure ||. It is described by two di- 
mensionless parameters, the peak height Uq/UbT and the 
asymmetry parameter a. When a — 1/2, the potential 
is perfectly symmetrical, while a = or 1 corresponds to 
maximal asymmetry. In addition to U{x), the full poten- 
tial $ contains a term proportional to the driving force F, 
which figures in the dimensionless group Fa/k^T. Thus, 
to specify our potential fully, we require the three dimen- 
sionless parameters Uo/kBT, a, and Fa/k^T^ as well as 
Dq and the repeat distance a, which set a time and a 
length scale. More generally, we expect that any form 
of U{x) with only one peak per period will be roughly 
characterized by a peak height C/q (equal to the differ- 
ence between the miniinum and the maximum values of 
C/(x)), and an asymmetry a (defined as the distance be- 
tween a minimum in U{x) and the next maximum to the 
right, divided by a). While we have no a priori infor- 
mation about a, we have suggested that C/q should be 
of order several fceT, and have argued Fa/k^T sa 5 for 
KBBD's experiments. 

Because the time required to diffuse over a barrier de- 
pends exponentially on the barrier height, small differ- 
ences in Uq can lead to significant changes in transloca- 
tion speed, consistent with KBBD's observations. Fur- 
ther, if U{x) is asymmetrical, forces F and —F will lead 
to different barrier heights, and thus to different mean 
drift speeds for the diffusing polymer. Figure |6| illustrates 
this point. It thus appears plausible that a polymer mov- 
ing through the pore with its 5' end first could travel at 
a very different speed from one going through the pore 3' 
end first. There is, however, one additional complication. 
As shown in figure M, three different "vector" quantities 
can be oriented relative to the membrane: the applied 
electric field, the a-hemolysin pore, and the DNA. Each 
can point towards the cis or the trans chamber. With, 
say, the pore orientation held fixed, there are thus four 
possible situations. The two that have been realized in 
the experiments of KBBD are related by a flip of the 
polymer, while transforming F ^-^ —F (or equivalently 
U{x) I— > U{—x) ) in our model amounts to changing the 
direction of the field. Thus, although a velocity difference 
of a factor of 3 can clearly be accounted for by reasonable 
choices of parameters, we cannot quantitatively address 
how the two peaks in KBBD's experiments are related. 
Once all four possible situations have been explored ex- 
perimentally, however, it should be possible, for example. 



to estimate the value of a by comparing data for the ap- 
propriate pairs of situations. 



Effective Mobility and Diffusion Coefficient 

We now turn to the task of calculating the parame- 
ters V and D that describe the behavior of Eq. g on long 
length scales. Several approaches are available; in this 
section, we will describe the results of an analysis based 
on ideas of Risken (Risken, 1984). Details of the calcu- 
lation, which relics on an eigenfunction expansion, are 
given in appendix B. In the most general case, v and D 
have fairly complicated forms, but relatively simple lim- 
iting cases capture most of the relevant behavior. For ex- 
ample, one finds (le Doussal and Vinokur, 1995; Scheidl, 
1995) 

i ^ J_ [°° dze^^^/'''''^ r '^^[U(x + z)-U(x)]/kuT ^ (^g^ 

V Do Jo Jo a 

from which a number of limiting behaviors can be ex- 
tracted. Several equivalent expressions for v, as well as a 
similar, but more involved, expression for D, can also be 
obtained. 

Figure plots the velocity v versus F for polymers 
traveling in two different directions in the same (asym- 
metric) potential. At typical values of F, differences in 
velocity between forwards and backwards motion of a 
factor of 3 or more are easily obtained. Likewise, the cal- 
culated velocities are much slower than they would have 
been in the absence of a potential. 

One can gain more quantitative insight into both of 
these observations by studying how v and D behave 
in various limiting cases. Relegating the derivations to 
appendix B, we next consider several such expressions. 
Three cases are particularly of interest: large and small 
driving force F, and large potential barriers C/q (the case 
of small Uq corresponds to the absence of a potential and 
was discussed earlier). For small F, v and D must satisfy 
an Einstein relation. Indeed, in this limit one finds. 



DoF 



1 



D = Do 



1 



1 



r(0) r(0) 
^1 -'2 



l + 0( 



Fa 



knT' 



(10) 



where 
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Thus, v/ D — F/kBT, as the fluctuation-dissipation the- 
orem requires, but the effective diffusion coefficient D is 
reduced from its bare value Do by a factor that grows 



exponentially with the characteristic height of the po- 
tential. Perhaps more surprising is the fact that a linear- 
response-like regime is also reached for sufficiently large 

F. As F ^ oo, 
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The physical content of this result is that when F is 
much larger than a typical force derived from U{x), 
^'{x) « —F, and contributions from U may be neglected 
entirely. In the opposite limit of large C/qj one might 
expect that the diffusion process can essentially be de- 
scribed as hopping from one potential minimum to the 
next. Approximate formulas based on the Kramers es- 
cape rate (van Kampen, 1992) should then apply. In fact, 
for large Uq one finds 
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As before, we select the origin of U{x) so that its maxi- 
mum and minimum in each period occur at points Xmax > 

■^inini Wltn 3^niax '^inin — OLd. 

We have already estimated from KBBD's data that 
Id = D/v w 40a. A striking feature of the asymptotic 
forms Eq. p^ through Eq. ^just obtained is that all three 
imply a much smaller value. As we noted when we intro- 
duced the parameter Id, the linear response results both 
yield Id = k^T / F; given our naive estimate Fa/k^T w 5, 
we find Id ~ a/5 ^ 40a. For Uq large enough that the 
"hopping" approximation of Eq. HJ applies, this order of 
magnitude is little changed even as i^ — > cxd. Indeed, in 
this limit Eq. O gives Id = a/2. It is of course possible 
that some particular form of U{x) with finite C/q and F 
might lead to a value of Id of order 40a. It seems more 
likely, however, that Id interpolates reasonably smoothly 
among its various limiting values. The inset to figure H il- 
lustrates this point for the sawtooth potential introduced 
earlier. Although v and D each separately can depend 
strongly on the shape of U{x), their ratio is far less sensi- 
tive. We are thus led to one of the central conclusions of 
this paper: While many aspects of KBBD's results can 
be qualitatively explained by a model of diffusion in a 
one-dimensional periodic potential, the observed width 
of their peaks is inconsistent with this model if one takes 
Fa Ki bk^T. 



DISCUSSION 

In the previous section, we argued that the peaks 
in KBBD's distribution of first passage times are 



much wider than is consistent with our minimal one- 
dimensional model. It is not difficult to suggest reasons 
why this might be the case. Perhaps the most obvious is 
that Fa/k-oT could differ significantly from 5. Not only 
would a decrease of a factor of 100 in F bring our pre- 
diction for Id into line with experimental observations, 
it would also explain the polymer's unexpectedly slow 
translocation speed. At least two effects might decrease 
F. First, unless the pore has infinite resistance, not all 
of the applied voltage drop V will be across the pore. Al- 
though the large {^ Mfi) resistance of the a-hemolysin 
channel makes it unlikely that this mechanism could di- 
minish F by orders of magnitude, it certainly leads to 
some decrease. Second, the fact that there is a nonzero 
ionic current flowing through the pore while the polymer 
is translocating means that the motion of the polymer 
itself need not satisfy detailed balance. That is, the er- 
ror rate, or ratio of the probabilities of moving forward 
one base to moving backwards one base, is no longer 
required to be equal to eyip[eV / k^T) . To use a some- 
what different language, as the counterions are forced 
through the pore by the electric field, they entrain some 
of the solvent along with them. This solvent fiow ex- 
erts an additional drag force on the polymer, and this 
drag contributes to the mean force F. As a result, the 
electrophoretic mobility of the polymer in the channel is 
not in general equal to its hydrodynamic mobility mul- 
tiplied by its charge. Appendix C presents simple esti- 
mates based on continuum mechanics that suggest that 
both of these effects are small. These estimates, how- 
ever, make a number of simplifications; indeed, even the 
validity of the continuum equations is not assured on the 
nanometer scale. Given the importance of a large value 
of Fa/k^T to any attempts to sequence polynucleotides 
using the a-hemolysin pore, it thus seems desirable to 
verify experimentally that it is indeed approximately 5. 
Although a smaller than expected driving force is cer- 
tainly one mechanism that would generate wider peaks, 
others exist that do not require a large error rate. In 
many ways, our most poorly justified assumption is that 
the motion of the polymer backbone through the pore 
is much slower than the relaxation of every other de- 
gree of freedom in the system, so we begin by consider- 
ing what might happen if this assumption were to break 
down. For example, the protonation state of the open 
a-hemolysin channel is known to fluctuate on a much 
slower time scale than the characteristic polymer time 
a/v ^ 1/xsec (Kasianowicz et al., 1995), and the en- 
ergy barrier to moving a base through the pore might 
change significantly when the protonation state changes. 
It is instructive to consider a naive extension of our one- 
dimensional model meant crudely to describe such a sit- 
uation. Suppose that the pore -I- polymer system can 
be in one of two states, state 1, in which the polymer 
backbone can diffuse freely, and state 2, in which the 
backbone is trapped and cannot move. Let there be a 
transition rate (per time) LOij from state i to state j. 
This situation bears some similarities to popular mod- 



els of motor proteins (Jiilichcr et al., 1997), but with 
the important difference that the ratio loyiIlo2\ need not 
violate detailed balance; a similar description has also 
recently been proposed for the one-dimensional motion 
of RNA polymerase along a polynucleotide (Jiilicher and 
Bruinsma, 1998). If Pi(x) is the probability that the sys- 
tem is in state i and that a length x of polymer has passed 
through the pore, the long time diffusion of the system 
is governed by equations of the form 
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In state 2, the motion is arrested, so both the velocity 
and the diffusion coefhcient vanish. Just as in the one- 
dimensional, periodic case (appendix B), this model leads 
to a spreading Gaussian wave packet, with velocity and 
diffusion coefficient determined by the behavior of the 
eigenvalues near zero. One finds a velocity 
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where Aw = 0^21 — tJi2 and lu — lui2+lu21- Thus, if wi 7^ 
and oj and Aw are chosen properly. Id = D/v can be made 
arbitrarily large. This is true even if Di/vi remains of 
order a. Thus, even in this simple example, broad peaks 
are possible as soon as one relaxes the constraint that the 
model only contain one degree of freedom. 

In addition to the possibility that there is more than 
one slow degree of freedom, many other factors could 
contribute to the wide blockage time distributions ob- 
served by KBBD. For example, the observed peaks could 
reflect not the distribution in passage times for poly- 
mers of a given length, but the length distribution of 
the polymers themselves. Although the polydispersity of 
the poly[U] samples used by KBBD was not well char- 
acterized, no qualitative differences were seen with a 
perfectly monodisperse sample of DNA (Dan Branton, 
Harvard University, personal communication), suggest- 
ing that polydispersity is not the culprit. Similarly, we 
calculated the free energy !F of the parts of the polynu- 
cleotide outside the pore using a model that applies to 
conventional polymers above the theta point. If JF took 
a different form, we might not be able to neglect it. In 
particular, significant asymmetries between the two sides 
of the membrane could result in a non-electrical contri- 
bution to the driving force (Park and Sung, 1998b; Di 
Marzio and MandeU, 1997; Carl, 1998). For example, 
if the polynucleotide adsorbs weakly on one side of the 
membrane, there would be a force towards the adsorbing 
side of order fk-oT/a, where / is the fraction of adsorbed 



monomers. A similar effect could be obtained by con- 
fining the polymer on only one side of the membrane. 
Indeed, for sequencing applications, it might be useful 
intentionally to introduce an asymmetry as a way to ma- 
nipulate the polymer's speed and error rate without af- 
fecting the ionic current. Similarly, the current can in 
principle be varied with little effect on the polymer by 
putting a high concentration of macroions (e.g. colloidal 
particles) on one side of the membrane. This would in- 
duce a concentration gradient in their counterions that 
tended to drive the ions across the membrane. 

Finally, we would like to touch on one other issue of 
particular relevance to efforts to sequence polynucleotides 
as they pass through the pore. All of our results up to 
this point have been strictly valid only for homopoly- 
mers. Since it is known that diffusion in random media 
can be qualitatively different from diffusion in ordered 
systems (Bouchaud and Georges, 1990), it is worth ask- 
ing whether we expect any important changes when the 
homopolymer is replaced by a random heteropolymer. As 
long as the assumptions leading to our model of diffusion 
in a one-dimensional, periodic medium hold, one can ar- 
gue that the effect of using heteropolymeric DNA would 
be to modify the potential U{x). Rather than having 
an identical form within each unit cell of length a, U{x) 
might take one of four different shapes, corresponding to 
four different bases. It is known that one-dimensional dif- 
fusion becomes anomalous when \U{x) — U{Q)Y ^ 00 as 
x -^ 00, where the overbar indicates an average with re- 
spect to the random distribution of bases (Bouchaud and 
Georges, 1990; le Doussal and Vinokur, 1995; Scheidl, 
1995). In biological DNA sequences, it is believed that 
successive bases are either uncorrelated or have correla- 
tions that decay algebraically with distance (Herzel et al., 
1998; Stanley et al., 1992). In either case, [U{x) - U{Q)Y 
remains bounded for large x. Thus, within the sim- 
ple one-dimensional model, no qualitatively new behav- 
ior would result from replacing homopolymers with het- 
eropolymers. This is what one would expect to observe 
for short polynucleotides like those used by KBBD. It is 
worth mentioning, however, that this conclusion is sensi- 
tive to other effects. For example, the electrophoretic mo- 
bility of the polymer in the pore presumably has some se- 
quence dependence; this would lead to an effective short- 
range correlated random force on the polymer. Such 
a random force, however small, would in principle re- 
sult in anomalous diffusion on sufficiently long length 
scales (Bouchaud and Georges, 1990; Fisher et al., 1998). 



CONCLUSION 

The central idea of this paper was first presented in the 
introduction: In the experiments of KBBD, and likely 
in other examples of the translocation of biopolymers, 
the channel through which the polymer passes cannot 
be viewed simply as a set of hard, homogeneous walls. 



Rather, more specific interactions between the polymer 
and the channel must be taken into account. Indeed, we 
have argued that there is a regime in which polymer-pore 
interactions dominate, allowing a quasi-one-dimensional 
description of the translocation process. One immediate 
consequence of this observation is that on long enough 
length scales, the transport of the polymer through the 
pore is governed by a simple phenomenological equa- 
tion. Starting from this equation, we have derived several 
predictions about the polymer's distribution of passage 
times. For example, we have shown that the polymer's 
mean translocation time depends linearly on its length 
only for an extremely long polymer. It is perhaps worth 
reemphasizing that none of these results depend on any 
particular microscopic model of the pore. In contrast, 
several important qualitative observations of KBBD can 
be understood in terms of a more microscopic picture in- 
volving a "tilted washboard" potential. The tilted wash- 
board model also lead us to point out that the distri- 
bution of passage times in KBBD's experiments was far 
broader than one might expect from simple estimates, 
or indeed from any model with only one degree of free- 
dom. We have suggested several ways that this discrep- 
ancy might arise. Some, such as a serious mis-estimate 
of the mean force F on the polymer, imply an error rate 
in DNA sequencing of almost 50%. Several others, how- 
ever, do not require revision of the estimated error rate. 
Most notable among these are polydispersity in the poly- 
mer lengths and a strong coupling between the polymer 
and another degree of freedom in the pore with slow dy- 
namics. Determining which mechanism is at work in the 
experiments of KBBD remains an important experimen- 
tal and theoretical challenge. 

Our conclusions suggest several experimental avenues 
that might be explored in KBBD's or some analogous 
system. Most obvious would be to try to measure the 
error rate (or equivalently, the driving force F) . At least 
a rough estimate of Fa/k^T might be obtained in several 
ways. Obviously, any experiment in which it is possible 
to detect the passage of a particular nucleotide through 
the pore gives one direct access to the error rate. Al- 
ternatively, at small enough applied voltage V, it should 
be possible to observe a linear response regime. Devia- 
tions from linear behavior would then be observed at a 
value of V such that Fa/kBT ~ 0(1). With linear re- 
sponse data from all four possible relative orientations, it 
should also be possible to observe that dv/dV is the same 
for pairs of peaks. More ambitiously, if one could exert a 
non-electrical force on the polymer strong enough that its 
mean velocity through the pore fell to zero, this would 
give a direct measurement of F. Such an experiment 
might be accomplished with modern micromanipulation 
techniques. 

Other experiments of interest might test the existence 
of a quasi-one-dimensional regime. With enough data 
on the length dependence of imax, for example, it should 
be possible to observe the predicted deviation from the 
simple guess imax fx L. Further, if this data could be 



extended to sufficiently long polymers, deviations from 
the curve of figure |j would provide information on the 
crossover to a regime in which the dynamics of the poly- 
mer outside the pore are slower than those inside the 
pore. Once the basic theory has been verified, a number 
of different directions remain open. For example, a study 
of the fluctuations of the ionic current with a polymer in 
the pore could provide evidence about whether there are 
important slow degrees of freedom other than the poly- 
mer backbone itself. Adding a time varying component 
to the applied voltage in the experiments of KBBD might 
provide advantages in sequencing applications; the char- 
acteristic frequency for motion from one base to the next 
appears to be a relatively low 10^ Hertz. The behavior of 
polymers in very narrow channels is a rich subject that 
has only begun to be investigated. 



APPENDIX A: DETAILS OF CALCULATION OF 
DISTRIBUTION OF BLOCKAGE TIMES 



We wish to solve the equation 
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where t' = vt, for P{x,t') subject to the boundary 
conditions that P vanish at a; = 0, L and the ini- 
tial condition P{x,t' = 0) = 6{x — xq). The right 
and left eigenfunctions of Idd^/dx^ — d/dx are respec- 
tively exTp{x/2ld) sin{knx) and exp(— a;/2?(i) sin(fc„x), so 
we have 

2 °° 
P{x,t') = - ^ g-A„t'g-^o/2i<i sin(fc„2;o)e^/2''* sin(fc„a;) , 

(A2) 

with kn = WK I L and A„ = Z^fc^ -I- 1/4/^. After taking a 
derivative with respect to x, setting x = L, and simplify- 
ing, one finds that the probability that the polymer will 
exit the channel at a; = L at time t' is 

X Y^ e~^"*'fc„ sin[fc„(a;o - L)] . (A3) 

n 

Note that for large L this is a very slowly convergent 
series. It is thus convenient to rewrite it using the 
Poisson resummation formula. This formula states that 
for any function /, Y.n=-^I{n) = Y.Z=-oo fi'^-^rn), 

where / is the Fourier transform of / defined by /(g) = 
J dx f {x)e^^^ . Rewriting the sum in terms of the Fourier 
transform of the summand moves L from the denomina- 
tor to the numerator of the exponential, yielding, after a 
little algebra. 
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With if written in this form, the first term is exponen- 
tially larger than all subsequent terms as L becomes 
large, facilitating the analysis of limiting cases. With 
the help of a few definite integrals a number of results 
can be obtained exactly. In particular, if one defines 
/(a) = /o°°d2//2/"exph(y- 1)2/472/], then 7(3/2) = 
7(1/2) = 2y^7r7. The first equality can be proven with 
the substitution z = 1/y, the second with the substitu- 
tion u = {y — l)/y/y; I {a) for other values of a may 
be obtained by succesive integrations by parts. Using 
these identities, one can show, for example, that the to- 
tal probability that the polymer will exit from the trans 
side {x = L) if it started at xo is 
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One can similarly obtain exact expressions for {t') and 
for higher moments. 

Thus far, we have allowed the polymer's starting point 
3^0 to be arbitrary. Since the polymer always starts en- 
tirely on the cis side of the channel, the case of interest 
to us is a;o -^ 0. Eq. (A5) shows that the probability 
that the polymer passes through the pore vanishes in 
this limit. This conclusion is, however, a pathology of 
our model. We can still obtain a meaningful conditional 
distribution of passage times (that is, a distribution of 
passage times for those polymers that do leave at a; = i) 
by normalizing ip by the total probability of passage. In 
the limit xq — > 0, one obtains 
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All the terms but the first are subdominant as L ^ cxd, 
and when they are dropped, we obtain (Q). Note that not 
only is Eq. the correct asymptotic form for large L, but 
also that all subsequent terms describe peaks centered at 
increasingly larger values of t' . Thus, even when these 
terms significantly modify the behavior of ip{t') as t' — > 
oo, they can have a very small effect in the vicinity of 



APPENDIX B: DETAILS OF CALCULATION OF 
MOBILITY AND DIFFUSION COEFFICIENT 

Exact Expressions 

In this appendix, we will derive the expression of Eq. O 
for the mean drift velocity w of a particle in a periodic 
potential and an analogous expression for the effective 
diffusion coefficient D. Since the linear operator £ (de- 
fined by Eq. H) is periodic, it must have eigenfunctions of 
the Block form ip^{k,x) = e^''''^u^{k,x), where |fc| < 7r/a 
and u^{k, x) is periodic with period a. The eigenfunc- 
tion ipni^^ ^) is defined by Cip^{k, x) = —Xn{k)i{}^{k, x). 
Since £ is a non-hcrmitian operator, right and left eigen- 
functions are not equal, so we distinguish between them 
with superscripts R and L. Likewise, the eigenvalues 
A„(fc) are not in general real. The eigenfunctions are 
labelled by a band index n and a wavevector k in the 
first Brillouin zone. If the polymer starts at a; = xq at 
t = 0, then P{x,t) may be expressed as an eigenfunction 
expansion: 

Pix,t) - V / dfce''=("-"'')ML(fc,a;o)*u^(fc,x)e-^"('=)* . 

(Bl) 

Because of the exponential decay e"'^"'-'^^*, the smallest 
values of ReA„(fc) determine the behavior of P at long 
times. One can prove that the lowest value occurs at 
Ao(fc = 0) = 0. Performing a saddle point injtagration 
about this point, one thus finds that as i ^ oojHl 

P{x,t) ~ ^ ^~{x~xo-vtf/{4Dt)^,R(,^ = 0,x) , (B2) 

v47r7)t 



where 



dAo 
dk 



and D 



d^Xo 



fe=0 



dfc2 



(B3) 



fc=0 



At long times P{x, t) is thus a spreading Gaussian, mod- 
ulated by a periodic function v^{k = 0, a;) that gives de- 
tailed structure on the scale of a. One can hence reason- 
ably interpret the constants v and D in the expression for 
the Gaussian envelope as the same constants that appear 
in the macroscopic Eq. g. 

In light of these expressions, we obviously want to 
study the behavior of Aq in the vicinity of /c = 0. To 
do this, it is convenient to rephrase the eigenvalue con- 
dition Ctp^{k,x) = \n{k)tp^{k,x) in terms of u^{k,x) 
as 



^^Note that this result is valid only for values of x such that 
the difference x — xq — vt <^ 0{\/t) for large t, so that the 
quantity in the exponential has a saddle point exactly at fc = 
instead of somewhere in the upper complex k half-plane. 
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Dn 



d 



+ ik 



dx 

= -Xn{k)u^{k,x) 






C + ik\2— 
ox 



d $'(x) 



A:rT 



-k' 



v^{k,x) 
Un{k,x) 



(B4) 



If we view the fc-dependent part of the operator on the left 
as a small perturbation on £, then finding the derivatives 
of Ao(fc) at fc = is formally the same as-,a problem in 
quantum-mechanical perturbation theory.O As usual, we 
pose the expansions 



\o{k) 



dk 



fc=0 



UQ{k, x) = Uq(0, x) + kuQ ' (x) + 
Ug (fc, x) — Uq (0, x) + fcup' (x) + ■ 



fc=0 



(B5) 



The "ground state" u^(k — 0,x) of the unperturbed 
problem can be obtained exactly by integrating Cu^{k = 
0,a;) = (Risken, 1984): 

u^{k = 0, x) = iVe-*(-)/'=BT _ ^ /■" dx' e*(-')/'=BT 

Do Jo 

(B6) 

The two constants of integration N and S are determined 
by normalization and by requiring that u^{k = 0,x) be 
periodic. Physically, S is the (constant) probability cur- 
rent density in the stationary state u^{k — Q,x). Since 



C 



Do 



_a_ 

dx 






^, Uo(fc = 0,2:) is clearly a con- 
stant, which we choose to be 1. The first correction to 
Xo{k = 0) is then just the "expectation value" of the 
perturbation in the "ground state" : 



dX 

dk 



fc=0 



dx{2 — 
ox 



<i>'{x) 
k^T 



u^${k = 0,a;) 



(B7) 



After a bit of algebra, one obtains 

n _ p-Fa/k^T\ 



V = Sa = Dna ^^ ;, r,. , 

a2[/i/2-(l-e-^"/'=-^)/3] 


(B8) 


where 




Jo a 


(B9) 


/2^ r^e-*(-)/'^--^,and 
Jo a 


(BIO) 


Jo a Jo a 


(Bll) 



"'"'"There is a slight difference in that in our case the pertur- 
bation to the original operator C has terms of order k as "well 
as of order k. 



Note that the integrals I\ and I^ defined in Eq. ITl| 
are just /i and I2 evaluated with F = Q. Finally, v may 
be put in the form of Eq. g by everywhere rescaling a 
as a I— !■ na and letting 71 ^ 00 (le Doussal and Vinokur, 
1995; Scheidl, 1995). This substitution is valid because 
we could have originally thought of U{x) as having period 
na instead of a for any positive integer n. 

To find the 0{k'^) correction to Ao(fc = 0) requires 
knowing the 0{k) correction to the ground state Ug ' {x). 
This is typically expressed as a sum of eigenfunctions 
of the unperturbed problem, which are assumed known. 
Instead, we will invert C by direct integration to find an 
analytic expressio n fo r Uq ' {x) . After substituting the 
expansions of Eq. B5 into the eigenvalue equation and 
equating terms of order k, we find 



Cuf^{x) 



ivuQ{k = 0, x) 



-iDn 



dx 



<^'{x) 
kBT 



u^{k^O,x). (B12) 



Everything on the right hand side of the equation is 
known, so Uq ' (x) can be found by integrating twice. 
The constants of integration are determined by demand- 



R,l 



ing that Uq ' (a;) be periodic and that its inner product 



R,l/ 



with Ug (fc = 0, x) = 1 vanish. Once Ug ' (x) is known, we 
can equate terms of order k^ to learn that 



D-Do = J ul;{k = 0,xyi(2-^ 



9 , *'(^)~\ R,l. X 

Ug (X) 



knT 



(B13) 



A certain amount of additional algebra finally leads to 
an expression for D: 



D _. , r hJi - (1 - e~P-/^^T),h 

Do hh-{l-e-P-/^-'^)h ' 

where the li are the same as before, 

'"' dx r dx' 



(B14) 



Ji = / — / e 

lo a. Jq a 



■S>{x)/kBT 



/(^'), 



(B15) 



J2 = / — [ — f '^ f,-H^)/kiiT ^<i>(x')/kBT f(^ii 

a ./n a ,/n a 















(B16) 



.3.r^r^/M, 



(B17) 



and 



5*0 



duf{k ^0,x) S 



/(^) = ^<(fc - 0,x) + --v.-^-,-, _ __ ^ (B^8) 



dx 



Do 



Just as with the expression of Eq. B8 for the velocity, we 
can find equivalent formulas for D by replacing a by na 
and letting n —f 00. 
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Approximate Expressions 

Having obtained exact expressions for v and D, we will 
now outline how the limiting forms of Eq. |lO| through 
Eq. n4 are derived. We will focus on the expressions for 
V] the manipulations required to obtain the analogous 
asymptotic forms for D are very similar. 

The behavior of v for large and small F is most 
easily studied starting from Eq. |^. Define A{z) = 
L dx/ae'Kp[U{x + z)/k^T — U{x)/k-BT], and note that 
yl is a periodic function of z with period a. Expanding 
A{z) in a Fourier series and Laplace transforming term 
by term as demanded by Eq. y leads immediately to the 
formula of Eq. Mfor the small F behavior; in particular, 
the leading behavior is determined by the constant term 
in the series, L dz/aA{z) — l\^ I2 ■ Similary, we may 
find the behavior of v for large F by successive integra- 
tions by parts: Dq/v = A{0)/F + A'{0)/F^ + ■■■. 

The starting point for finding the large U behavior is 
similar. Rewrite v as 



Dna 



da;e"^^^'/''^^e^^/''='^ 



X / dye 



U{v)/k^T^-Fy/kBT 



(B19) 



First the inner, and then the outer integral can then be 
evaluated by Laplace's method as U becomes large. Note 
that since U{x) is periodic, we must sum a geometric se- 
ries over an infinite number of extrema to find the asymp- 
totic expression for the inner integral. Also, the location 
of the first maximum depends on the lower bound of in- 
tegration X. Once the two integrals have been evaluated, 
an expression equivalent to Eq. [l4| fo llows immediately. 
The only difference is that in Eq. |14|, we have chosen to 
write /} and /j instead of their large U forms. 



APPENDIX C: ESTIMATE OF REDUCTION IN 
THE DRIVING FORCE 

In the discussion section, we mentionned two factors 
that should reduce the driving force F on the polymer 
from the naive value Fa/k^T « 5. In this appendix, we 
present order-of-magnitude estimates of how large a cor- 
rection these effects cause. The estimates are based on 
the equations of continuum mechanics in a simplified ge- 
ometry and use bulk parameter values. They thus neglect 
a number of subtleties, most notably the presence of sig- 
nificant numbers of charged groups on the a-hemolysin 
pore itself. 

We begin by looking at the additional drag on the poly- 
electrolyte due to the flow of oppositely charged small 
ions through the pore. Consider a cylindrical polymer 
of radius r inside a cylindrical pore of radius R, and as- 
sume that the distance 5 = R — r between the polymer 
and the pore satisfies 5 "^ r and 5 < k~^, where k~^ is 



the Debye-Hiickel screening length. Both conditions hold 
in KBBD's experiments. Then neither the steady-state 
density n of ions nor the solvent velocity u parallel to 
the cylinder axis vary too strongly with radial distance. 
In particular, there is no Manning condensation. In the 
presence of an applied electric field E (of order the ap- 
plied voltage V divided by the length of the channel) , a 
body force term enE must be added to the Stokes equa- 
tion that describes very viscous flow. Here e is the elec- 
tron charge, and we have assumed that the only ions 
present are monovalent cations; the figure we obtain for 
the additional drag will thus be an upper bound on the 
drag possible with ions of both signs. We may use the 
component of the Stokes equation along the cylinder axis 
to estimate 



enE ~ rj 



V4_ 

,52 



(CI) 



where rj is the solvent viscosity. Since the Stokes equation 
is linear, we have written the axial solvent velocity as 
u = Up + u\, where Ui is driven by the electrical force 
on the ions, and Up by the motion of the polymer. The 
electrical current density J due to the ions is then roughly 



J ' 



e nfiE 



en{up + Uj) 



(C2) 



Here fi is the mobility of a single counterion. The first 
term describes motion of the counterions relative to the 
solvent, and the second the convection of the counteri- 
ons by solvent motion. The value of J is set by the fact 
that the total current / ~ 2TTrSJ is known to be of order 
10 pA in KBBD's experiments; Up must be of order the 
polymer's translocation speed. Thus, the pro blem is re- 
duced to solving the two "equations" CI and C2 in the 
two unknowns n and Up. By substituting reasonable pa- 
rameter value s, o ne finds that the term proportional to 
Up + Ui in Eq. C2 may be dropped. Then, the extra drag 



force per length on the polymer is 



drag 
length 






/ 

fie 



(C3) 



Substituting a typical value of an ionic mobility (in bulk 
solution), and multiplying by the channel's length Zch ~ 
50A, we obtain a drag force of order 10~^dyne due to 
the ion-driven solvent flow. This is to be compared to 
bk-oT/a ft! 3 X 10~^dyne. We conclude that the presence 
of the counterions in the channel does not signiflcantly 
reduce F. 

Another factor that we should take into account con- 
cerns the assumption that any applied voltage drop V 
falls entirely across the channel. In reality, some electric 
field must "leak" out of the pore, decreasing the force 
that drives the polymer. The simplest model in which 
one can consider this problem assumes that the bulk so- 
lution is neutral and has a constant conductivity cr, so 
that the electrical current density and the electric field 
are related by J = ctE. This picture of course does not 
hold within a distance k~^ of the membrane, but should 
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be reasonably accurate on longer length scales. Focus 
only on one side of the membrane, so that the pore acts 
as a current source injecting a current / at the bound- 
ary of an infinite half-space. Far way from the pore, the 
current density, and thus the electric field, should decay 
like 1/r^. In the near field, one expects this decay to be 
cut off at a distance of order the pore radius R. Knowing 
the current /, we can immediately obtain the electric po- 
tential, and thus the voltage drop Vc across the channel, 
in terms of the voltage V imposed at infinity. Our rough 
estimate gives (in cgs units) 



V-Vc 



27r crr^ 



dr ' 



TTcri? 



(C4) 



The factor of 2 in front of the integral arises because we 
must include the effect of a voltage decrease on both sides 
of the membrane. Substituting experimental values for 
/, (T, and R, we find that this calculation suggests that 
Vc differs from V by only a few percent. 
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FIG. 1. Histogram of number of observed blockade events 
versus the lifetime of the blockade, for 210 nucleotide poly[U]. 
The numbers 1 through 3 label the different peaks. From 
KBBD (courtesy of Dan Branton, Harvard University). In- 
set: Typical time series of the current versus time in the ex- 
periments of KBBD, showing a transient blockade due to the 
translocation of a polymer (courtesy of Dan Branton, Harvard 
University). 



FIG. 2. The distribution ■i/>(t) of passage times plotted 
versus t for L/ld ~ 5. Both quantities are appropriately 
non-dimonsionalized, t as vt/L and i/'(f) as L%l}{t)/v. The 
dashed curve is a Gaussian with the same mean and variance 
as ^{t). 



FIG. 3. Plot of the relative width St/tmsiic of the peak in the 
distribution of passage times, versus Id/L. This curve may be 
used to obtain the quick estimate Id ~ 40 nucleotides for the 
system studied by KBBD. The das hed cu rve gives the L -^ oo 
asymptotic behavior, 5f/fmax ~ y^2ld/L. We have chosen to 
put Id/L instead of L/ld along the ordinate to allow smooth 
contact with this large L behavior. 



FIG. 4. vtmax/L plotted versus L/ld. Note that vtmax/L 
varies signficantly over the range of L/ld relevant to the ex- 
periments of KBBD, and in particular that it does not reach 
its asymptotic value of unity until well outside the range of 
this plot. Inset: Plot of tmax (nondimensionalized by Id/v) 
versus L (nondimensionalized by Id). The dashed line gives 
the large L limiting form L/v, the solid line the exact value. 
Note that although tmax appears to the eye to depend lin- 
early on L over much of the range of the plot, it still differs 
significantly from L/v. 



FIG. 5. Sketch of the sawtooth "cartoon" potential dis- 
cussed in the text. The potential has period a, and aa is the 
distance from one minimum to the next maximum. The pa- 
rameter Uo gives the energy difference between minimum and 
maximum. 



FIG. 6. Sketch showing how asymmetry in the potential 
can lead to different speeds for forwards and backwards mo- 
tion. A bias is applied to the unperturbed potential (A) at 
top so that it has the same average gradient in the two bot- 
tom pictures. The potential at right (B), however, has been 
reflected through the vertical axis before the gradient is ap- 
plied. It thus has larger barriers to hopping from one min- 
imum to the next than the potential at left (C), leading to 
slower dynamics. 
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FIG. 7. The four possible relative orientations of polymer, 
pore, and applied electric field. In KBBD's experiments, the 
relative orientation of the pore and field is fixed and the ori- 
entation of the polymer is allowed to vary, corresponding to 
cases B and D. In our microscopic model, the pairs (A, B) 
and (C, D) are related by _F i— > —F. 



FIG. 8. Plot of the (nondimensionalized) average velocity 
V from equation y versus the driving force Fa/k^T\ v is calcu- 
lated using our microscopic model with a sawtooth potential. 
The parameter values are Uo/kBT — 10; a = 0.7 for the upper 
curve and a — 0.3 for the lower curve. The potentials for the 
two curves are thus related by U{x) i—> U{—x). Inset: The 
diffusive length Id versus the barrier height Uq of the sawtooth 
potential, for fixed driving force Fa — SksT and asymmetry 
a — 0.7. Note that over the entire range of Uo, Id i$ a. 
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